The dimension of loop-erased random walk in 3D 
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We measure the fractal dimension of loop-erased random walk (LERW) in 3 dimensions, and 
estimate that it is 1.62400 ± 0.00005. LERW is closely related to the uniform spanning tree and the 
abelian sandpile model. We simulated LERW on both the cubic and face-centered cubic lattices; 
the corrections to scaling are slightly smaller for the face-centered cubic lattice. 
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INTRODUCTION 



Loop-erased random walk (LERW) is the path ob- 
tained from a random walk by erasing loops as they are 
formed, and was introduced by Lawler. The paths con- 
necting points in a uniformly random spanning tree are 
distributed as loop-erased random walk [1-3]. Uniform 
spanning trees in turn are closely related to the abelian 
sandpile model of self-organized criticality [I| , and prop- 
erties of loop-erased random walk manifest themselves 
in the avalanches of these sandpiles. In particular, the 
fractal dimension z of LERW is related to the scaling 
behavior of topplings in the sandpile models 0, @] . 



dimension of 3D LERW reference 



In two dimensions, LERW has dimension 5/4 [1|, 
and its scaling limit (with suitable boundary conditions) 
is known to be SLE2 @. In dimensions 4 and higher, 
LERW has dimension 2, though with a cube- root log 
correction in 4D [Ioj |. In dimension 5 and higher, the 
scaling limit of LERW is Brownian motion. By contrast, 
relatively little is known about LERW in 3D. 

Kozma proved that 3D LERW (in certain domains) 
has a sca ling limit that is invariant under rotations and 
dilations [11| . The Hausdorff dimension of 3D LERW is 
not rigorously known to be well-defined, but 1 < z < 
5/3 [12j. There have been a number of estimates of this 
dimension z, sometimes expressed in terms of v — 1/z or 
2v. Guttmann and Bursill estimated 2v — 1.600 ± 0.006 
in d = 2 and 2v = 1.232 ± 0.008 in d = 3 0. Bradley 
and Windwer estimated 2v = 1.571 ± 0.006 in d = 2 and 
2v = 1.230 ± 0.003 in d = 3 Q. Anton predicted v = 
8/13 exactly Ag rawal and Dhar estimated z(d = 

3) = 1.6183 ± 0.0004 Fedorenko, Le Doussal, and 
Wiese gave an expansion for spatial dimension 4 — e: 



z(d = 4-e) = 2 



Evaluating this series at e = 1, they estimated z(d = 
3) = 1.614 ±0.011 17]. More recently, Grassberger es- 
timated z(d = 3) = 1.6236 ± 0.0004 
Agrawal and Dhar's earlier estimate, 
even more precise: 1.62400 ± 0.00005 
are summarized in Table |U 



18j . contradicting 



Our estimate is 
These estimates 



1 < z < 5/3 (rigorous) Lawler [121 

1.623 ± 0.01 Guttmann & Bursill [13] 

1.626 ± 0.004 Bradley & Windwer [14] 

13/8 (conjectured exact) Anton [15] 

1.6183 ± 0.0004 Agrawal & Dhar [16] 

1.614 ± 0.011 Fedorenko, Le Doussal, & Wiese [17] 

1.6236 ± 0.0004 Grassberger [18] 

1.62400 ± 0.00005 present work 



TABLE I: Estimates of the dimension 2 of loop-erased random 
walk in 3 dimensions. 



SIMULATION DESIGN 

Most earlier simulations (with the exception of 
Agrawal and Dhar's simulations [lj|) look at the length of 
the loop erasure of a random walk run for a large number 
of time steps. The starting point of a loop-erased ran- 
dom walk has different statistical properties than typical 
points on the LERW path. For example, the winding 
angle variance at the starting point of LERW is different 

3 



than at a typical point 



To measure the dimen- 



sion of LERW, we would like to measure the length of an 
LERW path without the atypical starting point. 

Agrawal and Dhar's simulations instead created loop- 
erased random loops, and measured their lengths. Dhar 
and Dhar had argued that adding an edge to a tree cre- 
ates a loop of size t with probability £~ 2 ^ z , and that 
the next step of an LERW produces a loop of size I with 
probability s» £- 1 ~ 2 / z [jlj . Agrawal and Dhar estimated 
the LERW dimension z by looking at the sizes of these 
erased loops. In their estimate, it was necessary to ig- 
nore both small loops (because of lattice effects) and large 
loops (which were influenced by the LERW stopping con- 
dition), effectively leaving fewer length scales with which 
to estimate the dimension. 

In our simulations, we ran a random walk on an 
L x L x L torus while erasing contractible loops, until 
a noncontractible loop-erased random loop was formed, 
and reported its length. This random variable is equidis- 
tributed to the following: Generate a uniformly random 
directed subgraph of the torus, where each vertex has 



2 



out-degree 1, conditioned on there being no contractible 
cycles. Every vertex leads to a cycle that winds around 
the torus. Pick a random vertex, find the cycle that it 
leads into, and report its length. By taking many such 
measurements for different L's, we estimated the dimen- 
sion of LERW. 

The geometry of the torus affects the length distribu- 
tion of the noncontractible loop, but for different L's the 
effect is the same. While Agrawal and Dhar had to ig- 
nore both small loops and large loops, we only need to 
ignore the small L's (because of lattice effects), so it is 
easier to see the asymptotic behavior. 



CHOICE OF LATTICE 

Kozma proved that LERW on any 3D lattice converges 
to the same scaling limit 



which is invariant under 



dilations and rotations, provided that random walk on 
the lattice converges to Brownian motion. To measure 
the dimension or other properties of LERW, in addition 
to using the standard cubic lattice, we also tried the face- 
centered cubic (FCC) lattice. The FCC lattice arises 
from the densest packing of spheres in 3D. Each site is 
adjacent to 12 nearest neighbors, as opposed to 6 nearest 
neighbors in the cubic lattice. In this sense the FCC 
lattice is closer to being isotropic than the cubic lattice, 
and since the scaling limit of LERW is isotropic, we might 
expect that LERW on the FCC lattice behaves more like 
the isotropic scaling limit for smaller values of L than 
LERW on the cubic lattice. 

The most convenient way to simulate on the face- 
centered cubic lattice is to use the same LxLxL grid as 
for the cubic lattice, but with extra edges, where x and 
y are connected if x — y is one of (±1,0,0), (0, ±1,0), 
(0, 0, ±1), (±1, T-l, 0), (0, ±1, Tl), or (±1, 0, T-l)- In ad- 
dition to changing the lattice, this also changes the ge- 
ometry of the torus, making it a skew torus. But the 
geometry of this skew torus is plausibly better than the 
geometry of the ordinary torus, since the girth divided 
by volume 1 / 3 is larger. 

Our simulations suggest that the face-centered cubic 
lattice gives slightly better results than the ordinary cu- 
bic lattice for systems with the same side length L. It 
would be interesting to see what effect the choice of lat- 
tice has on simulations in higher dimensions, where there 
are lattices that are much better (by some measures) than 
the hypercubic lattice. 



LOOP HOMOLOGY 

In order for the measured loop length to be related to 
the dimension of LERW, we would need to know that 
the loops do not wind around the torus too many times, 
and in particular that the number of windings does not 





homology 


L = 64 


L = 1024 


L = 16384 




1,0,0 


0.6242 


0.6214 


0.6212 


CD 


1,1,0 


0.2804 


0.2813 


0.2814 


c lattic 
c torus 


1,1,1 

2,1,0 
2,1,1 


0.0613 
0.0163 
0.0084 


0.0618 
0.0169 
0.0087 


0.0619 
0.0170 
0.0087 




2,0,0 


0.0077 


0.0080 


0.0080 


o CJ 


2,2,1 


0.0007 


0.0007 


0.0007 




2,2,0 


0.0005 


0.0006 


0.0006 




3,1,0 


0.0002 


0.0002 


0.0002 




3,1,1 


0.0001 


0.0001 


0.0001 




3,0,0 


0.0001 


0.0001 


0.0001 




1,0,0 


0.8395 


0.8377 


0.8375 


.2 9 


1, 1, -1 


0.0781 


0.0784 


0.0785 


FCC latt: 
skew tor 


1, 1,0 
2,0,0 

2, 1, -1 

2, 1, 


0.0674 
0.0067 
0.0061 
0.0012 


0.0683 
0.0070 
0.0063 
0.0013 


0.0684 
0.0070 
0.0063 
0.0013 




1, 1, 1 


0.0007 


0.0007 


0.0008 




2, 2, -1 


0.0001 


0.0001 


0.0001 




2, 2, -2 


0.0001 


0.0001 


0.0001 



TABLE II: Empirical probability distribution of the homology 
(up to symmetry) of the noncontractible loop. 



grow with the side length L (as it would in 5 and higher 
dimensions). Therefore we measured the homology of the 
noncontractible loop in addition to its length, and found 
that in fact it does not grow with L (see Table [TTJ) . 

The fact that the loops with high probability do not 
wind around the torus more than O(l) times can be de- 
duced as follows. A classical result says that two indepen- 
dent random walks started distance R apart and run R 2 
steps in Z 3 have a constant chance of intersecting. Lyons, 
Peres, and Schramm 22j showed that whenever two inde- 
pendent random walks intersect with constant probabil- 
ity, the second random walk intersects the loop-erasure of 
the first with constant probability. Every time the ran- 
dom walk winds twice more around the L x L x L torus, 
there is a constant chance that the second time around 
it intersects the loop-erasure of the first time around. So 
the number of windings is stochastically dominated by a 
geometric random variable, uniformly in L. 



NOTES ON THE SIMULATIONS 

We collected an enormous amount of data (using high 
performance computing clusters), e.g., 10 9 data points 
for L — 16384 for two different lattices, where each data 
point requires « L 2 random walk steps, for a total of 
« 5 x 10 17 random walk steps. For the results to be 
meaningful, we require a high-quality random number 
generator, and we used one based on the advanced en- 
cryption standard (AES-256), which has been found to 
have excellent statistical properties 23J, \2A j . 
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lattice 


2 8 ,...,2 14 


2 9 ,...,2 14 


2 10 ,... 


,2 14 


2 11 ,... 


,2 14 


2 12 ,... 


>2 u 


2 13 ,2 


14 


cubic 


1 . 6 2 3 9 3 ±^0©etnf 


1 . 6 2 3 9 6 d^DOflTTT 


1.62397 ± 


.00001 


1.62398 ± 


.00002 


1.62398 ± 


.00003 


1.62399 ± 


00006 


-er-rT712 ± .00006 


-8TTT687 ± .00008 


0.11677 ± 


.00011 


0.11668 ± 


.00015 


0.1167 ± 


.0003 


0.1166 ± 


0005 




p = 3 X 10~ 22 


p = 0.03 


p = 


4 


p = 


75 


p = 


57 


p undct'i 


ned 


FCC 


1.62394 +JWfflTT~ 


1.62396 ± .00001 


1.62397 ± 


.00001 


1.62398 ± 


.00002 


1.62399 ± 


.00003 


1.62400 ± 


00005 


-OrnSoSi ± .00006 


0.09062 ± .00007 


0.09054 ± 


.00010 


0.09048 ± 


.00015 


0.0904 ± 


.0002 


0.0902 ± 


0005 




p = 3 X 10 -13 


p = 0.08 


p = 


41 


p = 


43 


p = 


56 


p undct'i 


ned 



TABLE III: Estimates of the 3D LERW dimension from the data. Different sets of system size L were used; the fit in the first 
data column was for L £ {2 s , 2 9 , 2 10 , 2 11 , 2 12 , 2 13 , 2 14 }, while the last column used only L € {2 13 ,2 14 }. The fits shown here 
are least-squares fits of log(E[loop length]) to functions of the form z log L + a, where we used 10 9 data points for each system 
size L. (The parameter a depends upon the lattice, but z is the same for both lattices.) The error bars given are the 95% 
confidence intervals (±1.96 standard deviations) of the fitted parameters (z and a). For each such fit we did a x 2 test, and give 
the p-value of the % 2 statistic. 



Following Agrawal and Dhar [l6} , we used a hashtable 
of points visited by the loop-erased random walk to iden- 
tify newly created loops. The storage requirements are 
then order 2, 1 - 6240 - rather than order L 3 . However, un- 
like the hashtables in other simulations 16,ll8|, rather 



than use linked lists in the event that two different points 
on the LERW accidentally hash to the same entry in the 
hashtable, we used an "open address" hash table, since 
they have less data structure overhead [25| . Open ad- 
dress hash tables are not normally suitable when entries 
can be deleted (such as when loops are erased), but in 
the case of LERW simulations, we may delete the points 
in an erased loop in reverse chronological order, and in 
this case open addressing works with deletions. 



ESTIMATE OF DIMENSION 

Table IIIII summarizes some simple least-squares fits of 
the data to estimate the dimension, which suggest an ex- 
ponent of 1.6240. In Fig. [T]we show our simulation data 
for the length of the noncontractible loop in a histogram 
type format. It is evident from this figure that the expo- 
nent 1.6240 is at least approximately correct. 

To estimate the exponent to four decimals from a plot, 
we need another way to present the data. Therefore we 
let Q q {L) denote the q th quantile for loop length on a 
system of size L. For example, (5o.s(16384) is the em- 
pirical median loop length for the torus of order 16384. 
For any q, we would expect Q q (L) to take the form 
a(q)L z for large L. More precisely, we would expect 
there to be a correction term, most likely of the form 

Q q (L) — a(q)L z + b(q)L v H , where y might be related 

to the exponent for a close encounter of the LERW path 
with itself. Then 



log(Q g (Li)/Q 9 (La)) 



b(q)/a(q) 



y-z 



tv- 

L 2 



") + ■ 



log(Li/L 2 ) \og{Li/L 2 ) 

In Fig. [21 we plot these ratios of empirical quantiles to es- 
timate the dimension. It appears that b(q) is negative for 
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(a) PDF of loop length (cubic lattice, cubic torus) 
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(b) PDF of loop length (face-centered cubic lattice, skew torus) 

FIG. 1: Probability density function for the length of the 
loop-erased random loop on the L x L x L cubic lattice (top) 
and face-centered cubic lattice (bottom). For both lattices, 
data is shown for L = 2 9 , 2 10 , 2 11 , 2 12 , 2 13 , 2 14 . For each L we 
effectively show a histogram of the loop lengths; the intervals 
are very short, but nonetheless contain many data points, so 
the histograms appear to be curves. Here the loop lengths 
have been scaled down by a factor of L 1 ' 6240 to compare the 
data for different system sizes. At this scale, the curves for 
the six different L's are indistinguishable, and appear to be 
one curve. The differences between the cubic and FCC curves 
arise because the FCC simulations were done on a skew torus. 
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(a) Estimate of dimension from ratios of quantiles (cubic lattice) 
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(b) Estimate of dimension from ratios of quantiles (FCC lattice) 

FIG. 2: (Color online) Here we plot the quantile ratios 
log(Q 9 (L)/Q 9 (L/4))/log(4) versus q for L = 2 11 , 2 12 , 2 13 , 2 14 
for the cubic lattice (top) and face-centered cubic lattice (bot- 
tom). The saw-tooth pattern, which is evident for L — 2 11 
and still visible for L — 2 12 , arises because the loop length is 
integer-valued. The saw teeth are more pronounced for the 
cubic lattice because it is bipartite. The curves for larger L 
become progressively more fiat, and appear to be converging 
to a value in the range 1.62400 ± 0.00005. 



small q and positive for large q, taking the value near 
q = 0.7. For this q the first correction term is 0, so the 
quantile ratios more precisely give z. The quantile ratios 
(as a function of q) appear to converge to a horizontal 
line at a geometric rate, in agreement with the formula. 
Additional analysis of the data suggests that the correc- 
tion exponent y — z is approximately in the range —0.8 
to —0.85. Based on these plots, together with the least- 
squares fits, we judge that the dimension z is likely to be 
1.62400 ±0.00005. 



TAIL BEHAVIOR 

Also of interest are quantities such as the probabil- 
ity that a loop-erased path is unusually short. In two 
dimensions it was recently shown that the probability 
that a LERW path is shorter than A times its expected 



leng th decays exponentially fast in (or faster) 

[26j |. We estimated the corresponding tail behavior for 
3D LERW, and found that it decays exponentially fast 
in A~ Q for a approximately 0.58 ± 0.02. Perhaps the 
correct exponent a is the reciprocal of the dimension. 
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